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Abstract 

Described  herein  are  methods  designed  to  obtain  accurate  broadband  surface  wave  disper¬ 
sion  measurements  on  two  spatial/frequency  scales:  continent-wide  (Eurasia;  20  -  300  s)  and 
regional  (Central  Asia,  within  20®  of  KNET;  5  -  50  s).  These  methods  are  bcised  on  well  devel¬ 
oped  frequency-time  and  floating-filter  analyses,  and  are  included  within  procedures  that  utihze 
relational  parametric  and  waveform  database  structures  with  rapid  graphics,  which  allow  mea- 
smrements  to  be  made  quickly  on  relatively  large  volumes  of  data  from  a  variety  of  source  regions 
recorded  on  heterogeneous  networks.  These  methods  cire  cmrrently  being  applied  systematically  to 
GSN/CDSN/GEOSCOPE/MEDNET  data  and  Kyrghyz  Telemetered  Seismic  Network  (KNET) 
data  to  yield  Rayleigh  and  Love  wave  group  velocity,  phase  velocity,  polarization,  and  amphtude 
measurements. 

The  complexity  of  Eurasian  structure  imphes  that  the  crucial  problem  is  to  extract  the  desired 
signals,  related  to  nearly  directly  arriving  waves  that  can  be  interpreted  deterministically,  from 
the  essentially  stochastic  interfering  multipaths  and  coda.  Continent-wide,  200  events  have  been 
processed  yielding  dispersion  measurements  for  more  than  3000  paths  to  date  (Levshin  et  ai, 
1995).  Regionally,  data  from  KNET’s  continuous  channels  have  been  accumulated  for  future 
einalysis.  All  waveform  (including  ‘cleaned’  waveforms)  and  peirametric  data  (including  dispersion 
mejisurements  md  frequency-time  images)  are  stored  in  CSS  v.  3.0  (with  appropriate  extensions). 

In  the  near  future,  experiments  will  be  performed  to  automate  this  method,  which  currently 
depends  criticeilly  on  humein  interaction.  The  resulting  measurements  and  cleaned  waveforms  wiU 
be  used  as  follows.  (1)  Group  «ind  pheise  velocity  measurements  are  being  inverted  for  2D  maps  of 
group  (Levshin  et  ai,  1995)  and  phase  velocities.  (2)  These  maps  will  be  inverted  for  hthospheric 
shear  velocity  structure  which  should  provide  significeintly  higher  resolution  tham  previous  models, 
especially  in  Centrail  Asia.  (3)  The  cleaned  waveforms  are  ideal  targets  for  surface  waveform 
fitting,  which  augments  and  complements  the  dispersion  measurements.  (4)  Higher  resolution 
phaise  velocity  maps  in  Central  Asia  will  be  used  to  develop  surface  wave  stacking/array  processing 
methods  for  regional  broaidband  eurays  located  in  geologically  complex  settings,  such  as  KNET. 
Such  methods  must  incorporate  differences  in  the  dispersion  characteristics  among  stations  in  the 
array. 
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Objective 


Because  approximately  two-thirds  of  Eurasia  is  surrounded  by  seismically  active  regions  (some 
of  which  extend  well  into  the  continent),  in  principle  Eurasia  provides  the  Earth’s  best  laboratory 
for  studying  continental  wave  propagation  and  for  the  use  of  seismic  wave  characteristics  to  map 
crustal  and  upper  mantle  structmres  on  a  variety  of  length  scales.  The  development  of  seismic 
networks  (e.g.,  GSN,  CDSN,  GEOSCOPE,  KNET)  in  the  Former  Soviet  Union,  China,  and  the 
coimtries  of  Southern  Asia  hold  out  the  promise  to  recdize  aspects  of  this  potential. 

There  are  two  meiin  objectives  of  the  studies  described  herein.  (1)  The  first  objective  is  to 
apply  well  estabhshed  methods  of  surface  wave  ancdysis  to  a  continent-wide  broadband  network 
(GSN,  GEOSCOPE,  CDSN,  MEDNET)  eind  to  a  broadband  regional  array  (KNET;  Kyrghyz 
Telemetered  Seismic  Network)  to  yield  Rayleigh  and  Love  wave  group  velocity,  pheise  velocity, 
polarization,  eind  ampHtude  measurements  on  two  spatial  scales  and  frequency  bands:  continent¬ 
wide  (Eurasia;  20  -  300  s)  and  regional  (Central  Asia,  within  20®  of  KNET;  5  -  50  s).  These 
measurements  will  be  used  in  the  futinre  to  invert  for  more  sharply  focused  lithospheric  models  on 
these  two  spatial  scales.  (2)  The  second  objective  is  to  use  the  estimated  regional  phase  velocity 
maps  to  develop  surface  wave  stacking/array  processing  methods  for  large-scale  broadbeind  arrays 
in  complex  geological  settings,  such  as  KNET,  to  attempt  to  reduce  the  magnitude  threshold 
for  the  analysis  of  intermediate  period  (5  -  30  sec)  sirrface  waves.  To  date,  prehminary  research 
results  have  been  achieved  only  for  the  first  objective. 

Preliminary  Research  Results 


Measurement  Procedures 

Problems  associated  with  the  estimation  of  accurate  surface  wave  chzuracteristics  (wave  velocities, 
amphtudes,  polarizations)  do  not  chcinge  in  nature  with  the  spatial  scEile  or  frequency  band  of 
interest,  although  they  do  change  in  magnitude.  The  most  significant  issues  concern  the  accrual 
of  high  quadity  data,  the  identification  and  extraction  of  unwanted  signals,  and  the  measmrement 
of  the  signals  of  interest. 

Data  qucJity  on  both  continenteil  (GSN/GEOSCOPE/CDSN/MEDNET)  and  regional  (KNET) 
scales  is  quite  good,  as  exemphfied  by  the  record  sections  shown  in  Figure  1.  The  main  problem  to 
be  faced  is  that  continents  are  complicated.  This  not  only  makes  interpretation  in  terms  of  struc¬ 
tural  models  difficult,  but  also  greatly  comphcates  measurements;  or  more  accurately  complicates 
the  identification  of  the  aspects  of  the  waveforms  on  which  measurements  are  to  be  apphed.  Omr 
aim,  then,  is  to  extract  the  signals  we  desire,  related  to  neeirly  directly  arriving  waves  that  can 
be  interpreted  deterministically,  from  the  potentitdly  interfering  multipaths  eind  coda  that  are 
essentially  stochastic  in  natmre.  Figmes  2  eind  3  demonstrate  the  analysis  procedure,  in  which 
imwanted  signals,  in  particular  surface  wave  coda,  overtones,  aind  body  waves  are  greatly  reduced 
in  the  filtered  seismogram  on  which  measurements  2ure  obtained. 

The  basic  characteristics  of  the  current  measurement  procedure  is  based  on  a  long  history  of 
development  of  surface  wave  cinalysis  (e.g.,  Dziewonski  et  al.  1969,  1972;  Levshin  et  ai,  1972, 
1989,  1992,  1994;  Cara,  1973;  Russell  et  al.,  1988).  The  recent  innovation  is  that  code  hais  been 
developed  which  allows  measurements  to  be  made  rapidly  on  relatively  large  volumes  of  data 
from  heterogeneous  networks  and  a  vziriety  of  source  regions.  The  innovations  have  required  the 
development  of  rational  parcimetric  and  waveform  database  structiares  (more  on  below)  and  the 
development  of  relatively  rapid  graphical  routines  for  human  interaction  with  the  data. 
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The  general  form  of  the  measurement  procedure  is  as  follows.  Group  velocity  -  period  diagrams 
for  the  vertical,  radial,  and  transverse  components  are  constructed.  An  analyst  manually  traces 
the  apparent  group  velocity  curve  for  the  Rayleigh  wave  (on  the  vertical  and  radial  components) 
and  the  Love  wave  (on  the  transverse  component).  Time- variable  filters  are  applied  around  the 
selected  c\n:ve  to  separate  the  desired  signal  from  the  ‘noise’.  This  results  in  fiiltered  group  velocity 
-  period  diagrams  for  which  contamination  from  interfering  signals  should  be  reduced.  Group 
velocity,  phase  velocity,  amphtude,  and  polarization  measurements  cire  automatically  obtained  on 
the  filtered  images.  Figures  2a  zind  3a  attempt  to  display  this  procedure. 

An  xmfortunate,  but  currently  still  necessary,  chairacteristic  of  this  procedure,  is  that  it  depends 
crucially  on  direct  human  interaction  with  potentially  large  voliunes  of  seismic  waveform  data. 
The  success  of  this  method  is  based  on  the  analyst  accurately  identifying  the  main  dispersion  ridge 
of  the  fundamental  modes,  separating  the  ‘direct  arrived’  from  surface  wave  coda  at  periods  below 
about  30  seconds,  inspecting  interpolation  near  spectral  holes,  and  truncating  the  measurements 
appropriately  at  long  periods  as  the  signals  weaken.  This  interaction  hmits  the  speed  with  which 
the  method  can  be  apphed,  and,  therefore,  the  volume  of  data  that  can  be  processed. 

To  date,  the  method  heis  been  applied  to  waveform  data  from  approximately  200  events  sur¬ 
rounding  Eurasia.  The  anadysis  of  significantly  Izirger  volumes  of  data  will  require  the  automation 
of  the  technique.  Attempts  at  automation  will  be  based  on  initial  continent- wide  group  and  phase 
velocity  maps  which  are  cmrently  rmder  development  (Levshin  et  al.,  1995).  This  will  increase  the 
volume  of  data  that  can  be  processed  and  should  result  in  improved  resolution  of  the  group/phase 
velocity  maps  and  the  resulting  structural  models. 

Measurements  on  a  Continental  Scale 

Due  to  the  high,  average  efficiency  of  surface  wave  propagation  across  Eurasia,  surface  wave 
measurements  can  be  made  at  periods  up  to  100  -  150  seconds  for  earthquakes  as  small  as  Mg  = 
5.0  that  propagate  across  the  entire  continent.  Of  course,  measurements  can  be  extended  to 
longer  periods  for  substantiadly  larger  events.  To  date,  event  volumes  have  been  accumulated 
for  GSN/GEOSCOPE/CDSN/MEDNET  data  from  approximately  300  events  occurring  between 
1988  md  1995.  Of  these  events,  shown  in  Figure  4a,  approximately  200  have  been  completely 
processed,  yielding  measurements  on  more  thm  3000  paths  (Levshin  et  al.,  1995). 

As  cin  exEunple,  group  velocity  measurements  for  a  single  station  (KEVO,  Finland)  for  one 
event  (Kuril  event,  10/9/94,  Mg  =  7.0)  are  shown  in  Figure  3b  for  the  Rayleigh  wave  (measured 
on  the  vertical  and  radial  components)  at  periods  between  about  20  and  300  seconds  and  for  the 
Love  wave  at  periods  between  about  30  eind  250  seconds.  Predictions  for  the  spherical  model 
PREM  (Dziewonski  and  Anderson,  1980)  are  shown  for  comparison. 

A  useful  by-product  of  these  analyses  are  ‘cleaned’  or  ‘filtered  seismograms’.  Figure  3b  shows 
a  compeuison  between  the  raw  and  filtered  seismogreuns  for  a  single  stationievent  pair.  Surface 
wave  coda,  overtones  and  body  waves  have  been  greatly  diminished  from  the  cleaned  seismograms, 
m2iking  them  an  ideal  tcirget  for  surface  wave  fitting  techniques  during  a  later  stage  of  this  research. 

Measurements  on  a  Regional  Scale 

On  a  regional  scale  at  shorter  periods,  smaller  events  {Mg  <  5.0)  cein  be  analyzed  similzurly. 
KNET,  situated  in  a  complex  tectonic  setting  in  Central  Asia  surromaded  to  the  East,  West  and 
South  by  significant  seismicity,  is  a  natural  site  to  focus  studies  of  regional  sccile  mezisruements. 
To  date,  approximately  200  events  (Figxue  4b)  with  Mg  >  4.0  within  approximately  20  degrees  of 
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KNET  have  been  extracted  from  KNET’s  continuous  channels.  Analyses  of  these  data  are  just 
beginning. 

Figure  3b  presents  an  example  the  Einalysis  of  these  data.  Seven  KNET  stations  were  operating 
during  the  passage  of  surface  waves  from  an  event  in  the  Qinghai  Province,  China  on  1/17/94 
(A  «  16  degrees,  Mg  =  4.8).  Rayleigh  and  Love  wave  group  velocity  measurements  are  shown  in 
Figure  3b.  Rayleigh  wave  measurements  are  quite  similar  across  the  array  above  about  20  seconds 
period  and  for  Love  waves  above  about  30  seconds  period  at  this  azimuth.  Variations  across  the 
array  at  shorter  periods  result  both  from  real  differences  along  the  various  wave  paths  necir  the 
network  and  also  from  Rayleigh  -  Love  interference,  which  can  be  significant  since  the  group 
velocities  of  the  two  wave  types  are  similar  in  this  period  range.  Cleaned  and  raw  waveforms  are 
presented  in  Figure  3c. 

Data  Base  Structure 

All  waveform  and  parametric  data  as  well  as  surface  wave  measurements  are  stored  in  the  CSS  v. 
3.0  relational  database  (Anderson  et  ai,  1990)  plus  extensions.  This  data  base  wiU  be  dehvered 
to  the  funding  agents  upon  completion  of  the  contracts.  The  standard  relations  (affiliation,  event, 
gregion,  instrument,  network,  origin,  sensor,  site,  sitechan,  sregion,  wfdisc)  are  augmented  with 
two  event  relations  modified  slightly  from  CSS  v.  2.8  (centryd,  moment)  and  three  extensions 
(disp,  ftdisc,  wfedit).  The  wfedit  relation  contains  information  about  the  time,  duration  cind  na¬ 
ture  of  waveform  problems  (e.g.,  clips,  gaps,  nonlineaiities,  interfering  events,  etc.).  The  disp  and 
ftdisc  relations  point  to  dispersion  measurements  zmd  group  velocity  -  period  images,  respectively. 
For  each  station:event  pair,  raw  and  filtered  group  velocity  images  are  output  and  pointed  to  by 
the  ftdisc  relation.  Dispersion  measurements  (group  velocity,  phase  velocity,  spectral  amphtude, 
polarization)  are  output  and  pointed  to  by  the  disp  relation.  Cleaned  or  filtered  waveforms  are 
output  and  pointed  to  by  a  cleaned  wfdisc  relation. 

Conclusions  and  Future  Plans 

The  abihty  to  make  relatively  rapid  and  accurate  surface  wave  measurements  (phase  and  group 
velocities,  amphtudes,  polarizations)  on  fciirly  leirge  volumes  of  data  has  allowed  us  to  begin  to 
process  systematiccilly  stuface  waves  propagating  across  Eurasia.  Dispersion  measurements  have 
been  made  for  approximately  200  events  on  a  continental  scale  recorded  at  GSN,  GEOSCOPE, 
CDSN,  and  MEDNET  stations.  Data  from  continuous  diannels  of  KNET  have  been  accumulated 
and  are  just  beginning  to  be  processed  in  support  of  a  regional  scale  study  of  Central  Asia.  These 
developing  data  sets  should  provide  significantly  higher  resolution  than  represented  in  previous 
models,  which  promises  more  accvuate  lithospheric  models. 

Futmre  developments  will  be  cilong  three  lines. 

•  To  increase  the  rate  of  surface  wave  processing,  measurement  methods  must  be  automated. 
This  will  be  based  on  the  use  of  preHminauy  group  and  phase  velocity  maps  resulting  from 
the  continental  scale  study.  See  Levshin  et  al.  (1995)  for  preliminaory  group  velocity  maps. 

•  Filtered  seismograms  will  be  used  for  surface  waveform  fitting,  which  augments  and  com¬ 
plements  current  measurement  procedures. 

•  Higher  resolution  phase  velocity  maps  in  Centred  Asia  will  be  used  to  develop  surface  wave 
stacking/array  processing  methods  for  regional  broadband  arrays  located  in  a  geologically 
complex  setting,  such  as  KNET.  Such  methods  must  incorporate  differences  in  the  dispersion 
characteristics  among  stations  in  the  stack. 
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(a)  GSN  Kuril  Event:  1 994 282 Ms=7.0 


YSSLHZ 
MAJO  LHZ 
PETLHZ 
MOJLHZ 
MA2LHZ 
HIALHZ 
YAK  LHZ 
ADKLHZ 
SSE  LHZ 
TATO  LHZ 
XANLHZ 
ENHLHZ 
LZHLHZ 
COL  LHZ 
CUZLHZ 
KMILHZ 
WMQ  LHZ 
LSA  LHZ 
CHTOLHZ 
BRVK  LbHZ 
ARULbHZ 
KEY  LHZ 
OBN  LbHZ 
ABKTLbHZ 
KONO  LHZ 
KIV  LbHZ 
GNILHZ 
ESK  LHZ 
ANTOLHZ 
TBTLHZ 
GRFO  LbHZ 
PABLHZ 


; 

a  ^ 

1 

A. _ A  . 

I 

- 

1 

1 

O8.‘O0K>0.000 

1994282 


06:1  OKK)^ 
1994282 


08:20KKI.OOO 

1994282 


08:30.‘00.000 

1994282 


06:40:00^)00 

1994282 


Oe:50H>OiXN> 

1994282 


09:00:00.000 

1994282 


09:10:00X00 

1994282 


(b)  KNET  Qinghai Event:  1994 27 Ms=4.8 
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Figure  1 .  (a)  Record  section  of  IRIS/GSN  vertical  waveforms  recorded  for  the  Kuril  event  of  10/9/95 
(Ms  =  7.0).  Waveforms  have  been  low  pass  filtered  with  a  lower  corner  of  20  seconds  period. 
Raypaths  are  inset,  (b)  Record  section  of  KNET  vertical  waveforms  recorded  for  the  Qinghai 
China  Event  of  1  /27/94  (Ms  =  4.8). 
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Figure  2.  (a)  Group  velocity  -  period  diagrams  for  the  vertical,  radial,  and  transverse  components  recorded 
at  the  GSN  station  at  Kevo,  Finland  for  an  event  in  the  Kurils  (10/9/94,  Ms  =  7.0,  A  =  58.5  degrees).  Est¬ 
imated  group  velocity  curves  are  shown,  (b)  Rayleigh  and  Love  wave  group  velocity  measurements  comp¬ 
ared  with  predictions  from  PREM.  (c)  Raw  (thin  solid)  and  filtered  (bold  dashed  lines)  waveforms  for  KEV. 
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Figure  3.  (a)  Group  velocity  -  period  diagrams  for  the  vertical,  radial,  and  transverse  components  recorded 
at  the  KNET  station  USP  for  an  event  in  the  Qinghai  Province  (1/27/94,  Ms  =  4.8,  A  =  16  degrees).  Estimated 


group  velocity  curves  are  shown,  (b)  Rayleigh  and  Love  wave  group  velocities  measured  across  the  array 
for  the  same  event,  (c)  Raw  and  filtered  waveforms  for  the  station  USP. 
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wave  magnitude.  by  body  wave  magnitude. 


